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Inductive measurements of third-harmonic voltage and critical current density in bulk 

superconductors 
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We propose an inductive method to measure critical current density Jc in bulk superconductors. 
In this method, an ac magnetic field is generated by a drive current /o flowing in a small coil mounted 
just above the flat surface of superconductors, and the third-harmonic voltage V3 induced in the coil 
is detected. We present theoretical calculation based on the critical state model for the ac response 
of bulk superconductors, and we show that the third-harmonic voltage detected in the inductive 
measurements is expressed as V3 = GsujIq / Jc, where uj/2-k is the frequency of the drive current, 
and G3 is a factor determined by the configuration of the coil. We measured the Iq-V^ curves of a 
melt-textured YBa2Cu307-5 bulk sample, and evaluated the Jc by using the theoretical results. 

PACS numbers: 74.25.Sv, 74.25.Qt, 74.25.Nf 



Large melt-textured RE-Ba-Cu-0 (where RE rep- 
resents rare-earth elements) superconductors can trap 
strong magnetic fields, and have been developed for high- 
field engineering applicationsi The trapped field is pro- 
portional to the critical current density Jc, and a nonde- 
structive and contactless method is desirable to measure 
Jc distribution in large-sized superconductors. 

Claassen et al? proposed a nondestructive and con- 
tactless method to measure Jc in superconducting films, 
and the scientific basis of this inductive technique was 
theoretically justifiedi^ In that method an ac magnetic 
field is generated by a sinusoidal drive current Iq cos ujt 
fiowing in a coil close to the upper surface of a supercon- 
ducting film, and the amplitude V3 of the third-harmonic 
voltage V3 cos(3u;t -I- iJa) induced in the coil is simultane- 
ously measured. For superconducting films of thickness 
df, V3 is detected when Iq > Ico, where Ico oc Jed/ is the 
threshold current at which the sheet current of the film 
reaches its critical value Jcdf and the ac magnetic field 
penetrates below the film.^ The Jc is determined from 
Ico^ and the current- voltage characteristics can be de- 
termined by measuring the frequency uj/2tt dependence 

Ofjc^ 

In this letter, we theoretically calculate the ac re- 
sponse of bulk superconductors based on the critical state 
modelfSii and derive the relationship between Iq and V3 
that are measured in the inductive method. Experi- 
mental setups and measured quantities in the inductive 
method for bulk superconductors proposed in this letter 
are basically the same as those in the previously devel- 
oped inductive method for films? Because the thickness 
of bulk superconductors is much larger than that for su- 
perconducting films, an ac magnetic field generated by a 
small coil placed above the bulk superconductors cannot 
penetrate below the lower surface of the superconductors. 
The harmonic voltages of bulk superconductors, how- 
ever, arise from irreversible vortex motion^ even when 
the ac magnetic field is shielded below the lower surface 
of the bulk superconductor. Our theoretical calculations 
described in this letter show that the Iq-Vz curves for 
bulk superconductors significantly differ from those for 
gjj^g^2j3j4j^ but are similar to those for superconducting 
rodsii^ 

First, we theoretically calculate the ac response of a 



bulk superconductor of thickness d situated aX —d < z < 
0, where the upper surface is at 2: = and the lower 
surface at z = —d. A single coil mounted just above the 
superconductor is used to generate an ac magnetic field 
and to detect the induced voltagei^ The bulk supercon- 
ductor is considered to extend infinitely in the xy plane, 
and thus the electromagnetic fields are not affected by 
the edges of the bulk superconductor sample. The edge 
effects are negligible if the distance between the sam- 
ple edge and the coil edge is larger than the radius of 
the coili^ We consider a coil of circular shape that pro- 
duces axially symmetric fields independent of azimuth 
9 in cylindrical coordinates {p, z) = {p,9,z). Current 
density J — J{p,z,t)6, sheet current K = K{p,t)6, 
and vector potential A = A{p, z, t)9 have only azimuthal 
components. Although the magnetic field H has both a 
radial component Hp parallel to the superconductor sur- 
face and a z component normal to the surface, we 
assume that Hp plays a crucial role in determining the 
ac response of the bulk superconductors. The shielding 
current Jg flowing near the surface prevents penetration 
of magnetic flux into the bulk superconductors, and Js 
mostly cancels out H-^ near the surface of the bulk su- 
perconductors, \Hp\ ^ \Hz\ for z « 0. 

A drive current /o cos ujt flowing in the coil produces 
an ac magnetic field, and its radial component at z = 
is given by 



Ho{p,t) ^ -/oFi(p) cosujt. 



(1) 



The coil-factor function Fi {p) is determined by the con- 
figuration of the coil. 
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where N is the number of windings, Sc — {R2 — Ri){Z2 ~ 
Zi) is the cross-sectional area of the coil in the pz plane 
at i?i < p < i?2 and Zi < z < Z2 where Zi > 0. 

Responding to the ac magnetic field Hq, shielding cur- 
rent density Js is induced near the surface of the bulk 
superconductor, — Aq < z < 0, where Aq is the flux- 
penetration depth that is smaller than the thickness d 
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of the bulk superconductor. The vector potential arising 
from Js is given by 



dz' ■ 



Jsip',z',t) 



{\P^P'\^ + \Z-Z'\2) 



where the integral region Sgc corresponds to the surface of 
the bulk superconductor in the xy plane, i.e., < p' < c» 
and < 9' < 2tt. The As is expressed as the multipole 
expansion, As — Agi + As2 + • • ■ for \z'\ ^ i?2, where 

^0 f ,2^1 Ks(p',t) 



A,i(p,z,t) ^ ^ I <fp 

^ Ssr. 



— (4) 

Ss. ■ (|p-p'|2+z2)l/2' 

As2{p,z,t) = — / dp -— . (5) 

The Ks{p,t) ~ j'^jy^ dzJs{p, z,t) is the sheet current, 

and Ms{p,t) = Ms{p,t)p = J°j^^^dz{zz) x Js{p,z,t) is 
the sheet magnetic moment. 

The magnetic flux linked in the coil of 1 wound 

by a thin wire is calculated as 



m = 



2ttN 



Z2 



pdp / dzA(p,z,t) 



LJo cos LUt + <S>sl {t) + <i>s2 (t) 



(6) 



where the first term in the right-hand side of Eq. © 
arises from the drive current, and Lc is the self- 
inductance of the coil. The remaining terms in Eq. JBJ , 
$si -|- + • • arise from the shielding current in the su- 
perconductor; <I>si is due to Kg and $^2 is due to Mg. The 
ratio of |<&s2|/|'l?si| is on the order of A0/-R2 for Z2 ^ i?2, 
whereas |$si| is on the order of LcIq. Because the ratio 
Ao/i?2 is typically on the order of 0.01 or less, we can 
neglect $s2 to obtain the Hnear response to Iq. The non- 
linear response of bulk superconductors, however, arises 
from ^s2, as we show below 

When the ac field does not reach the lower surface of a 
bulk superconductor (i.e., 2\Ho\ < Jed and Aq < d), the 
sheet current is given by Kg = 2Hq (where \Ks\ < Jed), 
the ac magnetic field at the surface of the supercon- 
ductor (z — 0) is also given by 2Ho,'^ and the mag- 
netic flux due to the sheet current is simply expressed as 
$si = —Lcilo coswt, where Ld is the mutual inductance 
between the coil and its image. ^ The $s2 is calculated 
as the line integral of Eq. © along the wire of the coil, 
and is given by ^s2{t) = -2npo dp pFi{p)Ms{p,t), 
where Fi{p) is given by Eq. jSJ. Induced voltage in the 
coil is given by V{t) = RcIq cos — d^/dt, which can be 
rewritten as 



V{t) — RcIq cos ut + uj{Lc 



2'Kpo / dppFi{p) 



Lcijlo sinust 
dMsip, t) 



dt 



(7) 



where Rc is the resistance of the coil. The amplitude of 
the last term of the right-hand side of Eq. lO is much 
smaller than the other terms, but the first and second 
terms do not contribute to the harmonic response. 

The shielding current Jg = JsiPi t)0 induced by the ac 
magnetic field of Eq. in the critical state modelSii is 



given as follows: Js{p,t) — ^Jc for —A^{p,t) < z < 0, 
and Js{p,t) = ±Jc for —Ao{p) < z < -~A^{p,t), where 
Jc is the critical current density. The signs '±' ('=f') 
represents ('-') for < < tt and ('+') 
for TT < < 27r»ii The fiux-penetration depths are 
given by Aq{p) = 2IoFi{p)/Jc and A±{p,t) = Ao(p)(l T 
cosa'i)/2. Then, the sheet magnetic moment Ms{p,t) = 

dz{'-z)Js{p, 2, t) is calculated as 



Ms{p,t) = ^JcAoipf (cos Ljt±^sin^ ujt 



[2llF,{pf/Jc_ 



4 ^ sin(2fc - l)ujt \ , , 



TT^ (2/c-3)(4fc2 - 1) 
Substitution of Eq. (jSJ into Eq. yields 



V{t) = FiC0s(wt-F^9i) 



Vn COS nuit. 



(9) 



n=2 



The harmonic voltages are given by = for even n 
2, 4, 6, ..piiandFn = bV^/^n^-^) for odd n = 3, 5, 7,., 
The third-harmonic voltage is 



G3 



Gzujil/Jc 
16^0 



dppF^{p), 



(10) 
(11) 



where G3 is the coil factor determined by the configura- 
tion of the coil. The fundamental voltage Vi cos{ijjt-\-'di) 
in Eq. ||5J) is affected by the impedance of the coil, Rc and 
u}{Lc — Lei). Equation (|10|l is similar to derived by 
Beani to determine averaged Jc over the surface of cylin- 
drical superconductors exposed to a uniform ac magnetic 
field. Note that Eq. (|10|l is applicable to evaluate local 
Jc distribution at the flat surface of large bulk supercon- 
ductors responding to a nonuniform ac magnetic field. 

Next, we demonstrate evaluation of Jc by using the 
theoretical result of Eq. (|10|l from our experimental data 
for a melt-textured YBa2Cu307-5 (YBCO) bulk su- 
perconductor sample. The YBCO sample was 12.1 x 
12.7mm2 and 0.9 mm thick. The third-harmonic voltage 
induced by the YBCO sample was measured by using 
the experimental setup described in Ref. i5||. Although 
our experimental setup was similar to that in Refs.|3 and 
'3, an additional cancel coil was needed to compensate 
for harmonic noise voltages.^ The coil used to generate 
the ac magnetic field and to detect the induced voltage 
had N = 400 turns and dimensions of inner diameter of 
2Ri — 1.0 mm, outer diameter of 2i?2 = 3.5 mm, and 
height of Z2 — Z\ = 1.0 mm. The distance between the 
coil and the surface of the sample was Zi = 0.25 mm. 
A dc magnetic field B^c was applied perpendicular to 
the flat surface of the sample (parallel to the c axis of 
YBCO). 

Figure n shows the square root of the measured third- 
1/2 

harmonic voltage as a function of the drive current 
Iq for the YBCO bulk superconductor sample. The flt- 
ted lines for the data of i?dc = 1 and 3T intersect the 
origin, and are consistent with Eq. (|10|l . We can evaluate 
Jc from the slope of the lines in Fig. ^ by using Eq. H1U|) . 
where the coil factor calculated by using Eqs. Q and 
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drive current, /q [mA] 
FIG. 1: Drive current Iq dependence of the square root 

1/2 

of the third-harmonic voltage Vg for a melt-textured 
YBa2Cu307_i sample at 77K and uj/2n — IkHz. A dc mag- 
netic field _Bdc was applied perpendicular to the wide surface 
of the sample. 

(fTT|l is G3 = 735r2-sec/m^ for the coil used in our mea- 
surements. The resulting values are 1.7 x 10^ and 
1.0 X 10^ A/m2 for B^c = 1 and 3T, respectively. 

The fitted line for the data of B^c = OT in Fig. ^ does 
not intersect the origin, suggesting the following empiri- 
cal equation for /0-V3: 

V3 = G3uiIo^Ici)VJc for/o>4i, (12) 

where Jc — 1.9 x lO^A/m^. The threshold current is 
Id = 11 mA, which corresponds to the ac magnetic field 
2fj,o\Ho\ ~ 2mT at the surface of the bulk supercon- 
ductor sample. The threshold current Id might be re- 
flected by the lower critical field or surface barrier field. 
The importance of such barrier fields has been realized 
in harmonic-voltage measurements^ 

A rough estimate of the electric field Eg induced at 



the surface of the bulk superconductor is given as the 
product of angular frequency u, ac magnetic induction 
2/xo|i?o|, and the field penetration depth Aq ^ 2\IIq\/Jc: 

E, - 2^ow|i/o|Ao - Wli/olV-^c (13) 

Both the Eg induced in the bulk superconductors and 
V3 induced in the coil are proportional to ujIq/ Jc, and 
thus, the ratio Va/i^s (~ 5 cm for the coil used in our 
measurements) is governed only by the configuration of 
the coil. The third-harmonic voltage of V3 0.05 mV, 
which is the typical order of Vis shown in Fig. ^ implies 
that the electric field induced at the surface of the bulk 
superconductor is on the order oi Eg 1 mV/m, which 
is much larger than that in the magnetization measure- 
mentsji^ Because Eg oc w, we can measure Jc-Eg curves 
(i.e., current- voltage characteristics) by measuring uj de- 
pendence of Jc, as demonstrated for filmsi^ 

In summary, we theoretically and experimentally inves- 
tigated the third-harmonic voltage V3 in response to an 
ac magnetic field generated by a sinusoidal drive current 
Iq in a small coil placed just above bulk superconduc- 
tors. The theoretical third-harmonic voltage induced in 
the coil is given by Eq. (|10|) . The experimental data of 
melt-textured YBCO bulk samples for B^c > 1 T are con- 
sistent with Eq. (|10|l . and the Jc values can thus be eval- 
uated from the experimental data. These results show 
that the distribution of Jc in large-sized bulk supercon- 
ductors can be determined by measuring I0-V3 curves 
with a scanning coil (or fixed multiple coils) and using 
Eq. H10|l . We can measure Jc only near the surface (i.e., 
Jc for — Ao < z < 0) of bulk superconductors, whereas 
information on Jc deep inside bulk superconductors can 
not be obtained by this method. 

This work was done as part of the Super- ACE project 
(a project for the R&D of fundamental technologies for 
superconducting ac power equipment) of the Ministry of 
Economy, Trade, and Industry of Japan. 
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